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We analyze equilibrium properties of coupled-doped cavities described by the Jaynes-Cummings- 
Hubbard Hamiltonian. In particular, we characterize the entanglement of the system in relation 
to the insulating-superfluid phase transition. We point out the existence of a crossover inside the 
superfluid phase of the system when the excitations change from polaritonic to purely photonic. 
Using an ensemble statistical approach for small systems and stochastic-mean-field theory for large 
systems we analyze static disorder of the characteristic parameters of the system and explore the 
ground state induced statistics. We report on a variety of glassy phases deriving from the hybrid 
statistics of the system. On-site strong disorder induces insulating behavior through two different 
mechanisms. For disorder in the light-matter detuning, low energy cavities dominate the statistics 
allowing the excitations to localize and bunch in such cavities. In the case of disorder in the light- 
matter coupling, sites with strong coupling between light and matter become very significant, which 
enhances the Mott-like insulating behavior. Inter-site (hopping) disorder induces fluidity and the 
dominant sites are strongly coupled to each other. 
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I. INTRODUCTION 

Quantum phase transitions are a remarkable zero 
temperature phenomenon driven by quantum fluctua- 
tions pp. Such transitions have been studied in many- 
body quantum systems where each quantum phase can 
be unambiguously defined. However, recent results show 
evidence that interesting aspects and important traces 
of the physics of novel quantum phase transitions may 
already be observed in the limit of very few interact- 
ing sites [2-6 . This is particularly clear in hybrid light- 
matter systems, such as coupled electromagnetic cavities 
doped with two level impurities, where a Mott insulating 
to superfluid crossover has been predicted for as few as 
six or seven sites [2 j. These systems feature a composite 
fermion-boson excitation in each site, hence the term hy- 
brid, and quantities like the variance in energy for each 
site have been used as markers for the transition between 
different phases [7]. However, entanglement, an unique 
quantum correlation with no classical analog which has 
been related to fundamental features of quantum phase 
transitions [8 , may be regarded as a more adequate or- 
der parameter [HE]. In this work, we study the system 
entanglement to show how such quantum correlations re- 
late to the behaviors the system may present. 

One possible Hamiltonian describing doped and cou- 
pled cavities is the Jaynes-Cummings-Hubbard (JCH) 
model [2j [10] which, in some limiting approximations, 
mimics the more typical and simpler Bose-Hubbard one. 
The similarities between both models have prompted the 
use of the latter as a basis for the analysis of quantum 
phase transitions in the former both for a large [11] and a 
very small number of sites [12 . However, the analogy to 
this simpler model ignores the internal structure of each 
site what prevents one from exploring the increased com- 
plexity of the JCH system. The implementation of such 



systems has been proposed in different quantum optical 
setups such as planar lattices of one mode cavities each 
containing one quantum dot [13] . photonic crystal mi- 
crocavities [13] . circuit quantum electrodynamics with a 
finite system approach [15 , and in trapped ions [16 . One 
of the greatest advantages of all these setups is the com- 
bination of highly controllable experimental conditions 
and the large effective size of each site that allows for 
the design of mesoscopic simulators of condensed matter 
systems. 

In many cases the JCH system is naturally disturbed 
by noise that usually takes the system out of equilibrium. 
However, even in equilibrium, disordered imperfections 
in the system preparation may induce transitions that 
drastically change quantum phases and their correlations. 
Disorder may manifest itself in very different and even op- 
posite effects. The lattice imperfections, that differ from 
site to site, may suppress quantum coherence inducing 
the spatial localization of quantum states destroying the 
system fluidity, which leads to compressible, despite non- 
fluid, glassy phases [T7] , However, disorder may induce 
fluidity under certain circumstances [6] [18] [19] . The hy- 
brid nature of the system also leads to interesting effects 
under the action of disorder, as we show in the following 
sections. 

In this paper, we show that the entanglement between 
different constituents of the JCH system can be used 
not only to characterize the already known phase transi- 
tions, also present in the Bose-Hubbard model, but also, 
and more importantly, to identify new behavior involv- 
ing the nature, either hybrid or bosonic, induced by the 
more complex JCH interaction. We address small and 
large quantum systems, extremes that present similari- 
ties and differences that are of great interest: while few 
sites are experimentally feasible in a controllable way, 
phase transitions are better defined in large samples. We 
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also analyze the entanglement and the disorder- induced 
effects of the JCH hamiltonian. For the analysis of the 
small system we enter deeply in the statistics of the en- 
sembles induced by disorder, since in principle a physi- 
cal observer could perform spectroscopic measurements 
of the system structure and obtain the disordered pure 
states (or at least quasi-pure) pertaining to the induced 
ensembles. For the analysis of the large system we re- 
sort to stochastic-mean- field-theory (SMFT), which was 
recently developed in [T9l [20] and allows us to study on 
site statistics. We show how the statistics of the system 
change under the various ways in which disorder may set 
in and also show the disorder-induced phase transitions. 

The analysis of the clean system is developed in section 
II with one subsection for the small limit and another for 
the large limit. The disordered small system is addressed 
in section III and the disordered large system is addressed 
in section V after a recollection of SMFT in section IV. 
Section VI concludes the paper. 



II. THE HYBRID SYSTEM: 
JAYNES-CUMMINGS-HUBBARD 
HAMILTONIAN 

The system studied here features a chain of sites each 
of which containing composite excitations, also known 
as polaritons, created by the interaction of a boson and 
a fermion. A typical experimental proposal for these 
systems is devised in resonant cavities, the bosons be- 
ing the photons that occupy the cavity mode and the 
fermions being two-level electronic transitions of the on- 
site dopants, as depicted in figure and described 
in [2j [21] . The Hamiltonian of these coupled doped cav- 
ities (with two level impurities) is the so called Jaynes- 
Cummings-Hubbard model and it is given by 




FIG. 1. Coupled cavities doped with two- level systems. The 
figure shows one possible realization of this system, in this 
case in photonic crystals, where cavities are defects in the 
periodic structure of the crystal. The two-level systems can 
be excitons in quantum dots or electronic levels of dopant 
atoms, for example. 

transition frequency of the dopant of the respective site, 
thus we define the ith site detuning = — uji. 

The polaritons are eigenstates of the intra-site 
(Jaynes-Cummings) hamiltonian and are given by 
\n+) = sm(Q n )\g)\n) + cos(# n )|e)|n — 1) and \n— ) = 
cos(0 n )\g)\n) - sin(0 n )|e)|ra - 1), with tan(2# n ) = 
—g^/n//\. The states \n) are photon number states and 
\g) and |e) are the ground and excited states of the dopant 
inside the cavity. Finally, the number of particles opera- 
tor in the ith site is given by Ni = a\di + crjai. 



A. Behavior of Small Sample Systems 



n 
(hj) 

with n being the number of sites, and Si being the intra 
site Jaynes-Cummings interaction between the dopant 
and the resonator 

Si = Uia\ai + Vio\<Ti + 9i{o\ai + o^at). (2) 

The ith site annihilation operators are a$ and <Ji for the 
bosonic and fermionic species, respectively. In Eq. ([!]), 
T(ij) describes the photon hoping, or tunneling, between 
nearest neighboring sites 

T (iJ) = ~ A (hj) i4 a j + a i a )\ ' ( 3 ) 

The coupling strength between the two level system and 
the cavity in the ith site is given by and the photon 
tunneling strength between nearest cavities is A^^+iy 
The photon frequency at the ith. site is uji and Vi is the 



Recent works show that the system described in the 
last session undergoes a Mott-superfluid phase transition 
when going from small hoping to large hoping or from 
negative detuning to positive detuning [2j[10]. In the first 
case the transition is induced because the hoping strength 
circumvents the photon blockage regime (non-linearity 
due to the dopant-cavity interaction) and in the second 
case the excitations are directly driven from mainly elec- 
tronic (electrons are not able to hop) to mainly photonic, 
hence the fluidity. This phase transition can be witnessed 
by single site properties. For example, when the system 
is isolated and the average occupation number per site 
(Ni) is one (same number of excitations and sites), the 
variance of Ni for any given site is a good order parame- 
ter [2j |3] : in the Mott phase each site has a single particle 
and there is no number fluctuation whereas in the super- 
fluid phase the onsite number variance is maximum. This 
analysis begins to fail when one takes into account in- 
teractions with the environment and spatial fluctuations 
that may not preserve the total number of particles in the 
system. For instance, dissipation introduces variances of 
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the occupation number in each site and var(A^) may over- 
estimate fluidity. Furthermore, although the measure of 
var(A^) hints at the type of excitation that dominates 
each phase it cannot reveal this fundamental property 
in detail because it does not fully distinguish between 
photons, electrons and polaritons. 

We proceed to show that the entanglement between 
different constituents not only reproduces previous re- 
sults but actually allows for the identification of a new 
crossover that the previous analysis did not reveal. 




FIG. 2. Results for a two site system. (Top-a) Entangle- 
ment between the sites. (Middle-b) In-site Entanglement. 
(Bottom-c) Atom- atom entanglement. 

In order to quantify the entanglement between the var- 
ious components of the system we choose the negativity 
measure [22 , which is very convenient to calculate. It 
should be reminded that null negativity does not neces- 
sarily imply null entanglement, in fact, null negativity 
means null or bound entanglement. However, any non 
zero value of negativity guarantees some form of distil- 
lable entanglement in the system what will prove to be 
enough for distinguishing the different quantum phases 
of the system as a whole. Given the quantum state of any 
two constituents A and B of the system, their negativity 
can be found by partially transposing their reduced den- 
sity matrix R = pj^ and then summing up the moduli 
of the negative eigenvalues of R. 

We begin with the case of small systems, where the di- 
agonalization of the Hamiltonian is computable, by look- 



ing at properties of the lowest energy state \G) with the 
constraint of having equal number of excitations and sites 
T~L\G) = E n \G). In other words, \G) is the lowest en- 
ergy eigenstate of the Hamiltonian having n = (Y™ A^), 
with N{ the number operator at the ith site. Con- 
sider, now, a cluster of two sites, that is the small- 
est possible such system. Even for this very basic unit 
cell, the entanglement between the sites clearly presents 
the signatures of Mott and superfluid phases that were 
found for much larger systems in previous works. In 
the Mott phase (with one polariton per site) there is 
no entanglement between sites with the Mott insulat- 
ing state being |G(MI)) = 1 1 — > 1 1 — > (for two sites). In 
the superfluid phase and when the excitations become 
mainly photonic the sites become entangled, with the 
superfluid state for two sites (described in [3]) given by 

|G(SF)> = \g)\g) [^|11> - §(|20> + |02»] . It should be 

kept in mind that for the finite system analysis there is no 
phase transition, only a smoother crossover, even though 
the phase transition terminology is commonly adopted. 
The behavior of the system (phase-like diagram), quanti- 
fied by the entanglement between different constituents, 
is depicted in figure Q , where we show the entanglement 
between sites, the in-site entanglement and the entangle- 
ment between atoms. 



The site-site entanglement shows the phase crossover 
as partially presented in [4] . When the site-site entangle- 
ment is negligible the system resembles a Mott-insulator 
and when the site-site entanglement is non-negligible the 
system presents superfluid-like behavior. Thus the site- 
site entanglement indicates the regimes in which the sys- 
tem is insulating and superfluid with small and large val- 
ues of entanglement, respectively [figure (|2ja)]. In order 
to quantify the polaritonic behavior we can look at the in- 
site entanglement that measures how correlated are the 
photonic field and the electronic transition in a given cav- 
ity (or a site) [figure (|2jb)]. In the Mott-like regime (small 
site-site entanglement) the in-site entanglement is signif- 
icant, indicating that the system presents polaritonic be- 
havior. Deep in the superfluid-like regime (large site-site 
entanglement) the in-site entanglement is small, indicat- 
ing predominant photonic behavior. However, during the 
crossover (as a function of either A or A) entanglement 
presents a non-monotonic behavior, with a region where 
it is maximum. Such non-monotonic increase, which is 
even more pronounced in the atom-atom entanglement 
[figure (|2jc)], suggests that as the system size increases 
and reaches the thermodynamic limit a phase transition 
should be verifiable, i.e. since at the point of the phase 
transition there are fluctuations over all length scales, 
more degrees of freedom interact with each other such 
that entanglement can exist between more degrees of free- 
dom. Furthermore, the regime in which in-site and site- 
site entanglement coexist corresponds to a polaritonic- 
superfluid, rather then just a photonic-superfluid. 
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Large Sample Systems and the 
polariton-photon crossover 



We can also obtain a wider view of the system phase 
diagram varying the number of polaritons in the system. 
In order to do that we can couple the system to a chem- 
ical reservoir of polaritonic particles with chemical po- 
tential ji (at zero temperature), such that the system is 
in equilibrium with this reservoir. The chemical poten- 
tial can be explicitly included in the system Hamiltonian 
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FIG. 3. (Top) The Mean-field parameter. (Middle) Entan- 
glement between the sites in the cluster. (Bottom) The in-site 
entanglement, that is, the entanglement between the dopant 
and the oscillator mode in mean field theory. All data with a 
four-sites cluster with the dimension of each oscillator trun- 
cated to 6 photons and zero detuning (A = 0). 



for the cluster becomes 

(cluster) 

+ + - \ aj \ 2 ], 



(i(cluster) ,j (environment)) 



(4) 



with (a) = a being the mean-field order parameter that 
has to be self-consistently determined by minimizing the 
ground state energy The phase diagram as a function of 
the chemical potential and the hoping frequency is shown 
in figure Q for large systems. 

Varying the chemical reservoir we can see the Mott 
lobes (each lobe corresponding to plateaus of different 
integer numbers of polaritons) in the infinite system [fig- 
ure ([3])]. The mean field parameter is null in the Mott 
phase and is positive in the superfluid phase. Only in 
this case we compute a site purity (one minus the purity 
more precisely) as the estimate of the entanglement be- 
tween such site and the rest of the chain (in this case, the 
cluster). Although this entanglement is not strictly zero 
in all of the Mott phase it still gives a fair account of the 
phase diagram and the lobe structure. We have consid- 
ered a four-sites cluster, which is a rather small cluster, 
even though it already requires a considerable computa- 
tional effort. Larger clusters would increase the precision 
of the site-cluster entanglement. The site-cluster entan- 
glement is maximum in the lobe borders (middle panel 
of fig. (|])), which, again, indicates strong polaritonic 
fluidity in the vicinity of the phase transition. 

Now, looking at the in-site entanglement (bottom of 
fig. ([3])), which can be regarded as the very essence of the 
polaritons, we can see the whole picture with the over- 
lay of the Mott- Superfluid and Hybrid-Boson crossover. 
The highest in-site entanglement is in the Mott-lobes and 
the lobe structure can also be defined by this quantity. 
Outside the lobes fluidity sets in, however the in-site en- 
tanglement is still very high indicating that the system 
has not yet undergone the Hybrid-Boson crossover de- 
spite having changed from insulating to superfluid. Far- 
ther away from the lobes and deeper into the fluid phase 
we finally observe the in-site entanglement vanishing in- 
dicating that the system finally turns bosonic. 



III. DISORDERED SMALL QUANTUM 
SYSTEMS 



For large (infinite) systems we adopt the mean-field 
approach, in which we treat a small cluster of sites in- 
teracting with a mean field, that is, a classical approx- 
imation of the rest of the chain (to which we refer to 
as an environment). This approach gives a factorable 
approximation of the non-factorable tunneling (or hop- 
ing) term by approximating the operators for their mean 
values plus a small fluctuation (in this case a quantum 
fluctuation) a = (a) + 5a. The mean field hamiltonian 



Every physical system presents imperfections (disor- 
der), that is, the system parameters may vary from site 
to site. There are many possible origins of disorder, for 
instance, imprecisions in the system manufacturing pro- 
cess, thermal fluctuations, and fluctuations induced by 
other uncontrollable electromagnetic sources in the sys- 
tem environment. One way to study the effect of disorder 
is to describe the parameters of each site as a stochastic 
variable & and the Hamiltonian becomes dependent on 



5 



the stochastic parameters Naturally, the system 

ground state becomes dependent on the values assumed 
by the system parameters \G) — >• and there 

emerges a new state, an average state 

P = J dp({ti})\G({ti}))(G({£i})\, (5) 

that contains the statistics of the effects induced by the 
static disorder, with being the distribution mea- 

sure of the disorder, such that it gives all moments of the 
site parameters = j dp • We choose to analyze 
only uncorrelated disorder such that the global measure is 
a product of local measures dp^i}) = Yl i with 

gaussian distributions P(£i) = ^J^ s exp{— }• The 

magnitude of disorder is then given by the distribution 
width or the mean square deviation S. 
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FIG. 4. The ensemble average entanglement E[p] between the 
two sites under disorder in the matter light detuning A (more 
specifically in the light frequency) on the phase diagram with 
5(A) = 10g. 

We can then characterize the average properties of the 
system given that it presents disorder. For instance we 
can look at the entanglement description of the phase 
diagram, only now we average the entanglement over 
the pure state ensemble generated by the different val- 
ues assumed by the system parameters. We remark 
that the pure state ensemble given in equation ([5| is 
a physically realizable ensemble [24], since in principle 
the experimentalist can perform spectroscopic measure- 
ments and obtain the values assumed by the parameters 
in that particular sample system and then prepare the 
system ground state. We can define the reduced states 
Pab({&}) = toE{\G({ti}))(G({Zi})\}, with the trace be- 
ing performed over the environment of A and B. For in- 
stance, if we are looking at the atom-atom entanglement 
then we trace out the field, so the field would be the 
environment in this case. Therefore, we can define the 
average entanglement between any constituents A and B 

E[Pab] = / dp({^})E[p AB ({^})}, (6) 



which is physically realizable since the ensemble of 
ground states also is [25j [26] . Notice that the aver- 
age entanglement of the ensemble is in general differ- 
ent than the entanglement of the average state with the 
usual hierarchy E[pab] > ^[pab], with the average state 
Pab = /dp({&})pAB ({&})• 
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FIG. 5. Probability distributions induced by disorder in 
the matter light detuning. (Top) Renormalized histogram 
P((Ni)) of the site occupation number as a function of disor- 
der. (Bottom) Renormalized histogram P(E) of the site-site 
entanglement as a function of disorder. The light hoping is 
A — g and average detuning (A) = 5g. 

In what follows in this section we consider only two 
sites of the jaynes-Cummings-Hubbard hamiltonian. 



A. Disorder in the matter- light detuning 

Now we can describe the effects induced by disorder in 
each of the parameters individually, detuning A, hoping 
A, and matter light coupling g. Let us begin by analyzing 
disorder only in the cavity- atom detuning, thus {&} = 
{A^} (see figures Q and (J5|). As can be seen in fig. Q 
the average entanglement between sites seems to decrease 
over the whole phase diagram in comparison with the 
clean case of fig. 

The decrease of site-site entanglement indicates that 
the excitations tend to localize through an Anderson-like 
mechanism. For instance, starting with the system in the 
superfluid phase as we increase the detuning disorder the 
distribution of the single site number occupation P((Ni)) 
is broadened and then it becomes a two peaked distribu- 
tion (see top panel of figure In the regime in which 
the distribution P((Ni}) presents two peaks the system 
is fully localized, such that one of the peaks corresponds 
to all excitations in cavity one and the other corresponds 
to zero excitations in the cavity. This extreme regime 
of localization can be regarded as a bosonic bunching: 
the large disorder in the cavity line width allows for re- 
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FIG. 6. Effects of static disorder in the light hoping A on 
the phase diagram. The ensemble average entanglement E[p] 
between the two sites. Disorder of S(A) = lOg. 



alizations in which the cavity has a very low frequency 
such that it is energetically favorable to fit more than one 
excitation in one cavity instead of distributing the exci- 
tations over the sites. For two sites and two excitations 
the state can be approximately given by |2)|g)|0)|g) or 
|0)|g)|2)|g) with the atoms in their ground states. 

The parameters in fig. (5) are such that the system is 
in the superfluid phase in the clean limit. In this case, as 
expected, the distribution of the site-site entanglement 
P(£'[SS]) shows a peak at the maximum value (bottom 
panel of the figure). However as disorder increases there 
emerges a second peak close to the minimum value corre- 
sponding to localized states. Thus the ensemble presents 
both superfluid and insulating states for intermediate val- 
ues of disorder. The presence of the two sorts of states 
can be regarded as a precursor of a glassy phase [5 , which 
we will show to be true with the large system analysis. 
Furthermore, as disorder increases even further the sys- 
tem becomes fully localized and the site-site entangle- 
ment distribution becomes single-peaked at very small 
values of entanglement. 



B. Disorder in the photon hopping 

Disorder in the photon hoping generates a different ef- 
fect (see figures (|6| and [7]). Carrying out the same anal- 
ysis as for disorder in the detuning, we see that the aver- 
age site-site entanglement increases in the region where 
a Mott phase exists in the clean limit, while it remains 
practically unaltered in the superfluid phase. Fluctua- 
tions in the hoping may actually induce a glassy fluid 
phase [6j[T8l[T9], that is, disorder allows realizations in 
which the hoping is stronger than the photon blockade 
and those realizations may prevail. The A-disorder may 
also suppresse the polaritonic behavior which can be seen 
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FIG. 7. Probability distributions induced by disorder in the 
light hoping. (Top) Renormalized histogram P((Z)) of the 
total atomic occupation number as a function of disorder. 
(Bottom) Renormalized histogram P(E) of the site-site en- 
tanglement as a function of disorder. The average light hoping 
is (A) = lg and detuning A = —2g. 

as a decrease in the atomic population. We can also 
look at the distribution for the total atomic excitation 
Z = "^ Ji (j\(J% as a function of disorder starting from the 
system at the Mott phase in the clean limit. The dis- 
tribution P((Z)) is very asymmetric and as disorder in- 
creases it concentrates at the extreme values assumed by 
(Z). One of the extremes corresponds to polaritonic su- 
perfluidity and the other to photonic superfluidity, with 
the latter prevailing in the limit of very large disorder. 
This can be corroborated by the entanglement distribu- 
tion P(£[SS]). 

C. Disorder in the matter-light coupling 

Finally we analyze disorder in the matter light cou- 
pling g (see figure (J8|). A first look at the ^-disordered 
average entanglement (not shown here since it resembles 
very closely the A-disorder case) suggests that disorder in 
the Jaynes-Cummings coupling also induces localization, 
that is, disorder allows a great number of meaningful real- 
izations in which the sites are almost unentangled and the 
excitations tend to bunch. However, the ^-disorder in- 
duced distribution of the site occupation number P((Ni)) 
(top panel of fig. Q) is very different than the one in- 
duced by A-disorder (which we showed that presents lo- 
calization). In the current case the P((Ni}) distribution 
shows three peaks, the extreme ones corresponding to 
bunching similar to the A-disorder case, and a middle 
one that corresponds to states in which the excitations 
are still equally distributed among the sites. Disorder 
in the matter-light coupling induces states with Mott- 
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like features in which the site-site entanglement vanishes 
(see middle panel of the figure). In fact, this is the 
meaning of the middle peak in P((Ni}): some sites un- 
dergo a superfluid-insulating transition through a Mott- 
like mechanism. Since the distribution of the atomic oc- 
cupation P((Z)) (bottom panel of fig. ^) is narrowly 
centered at an appreciable (although not extreme) value 
we may conclude that the system nature becomes mainly 
polaritonic, and thus presenting both Mott states (middle 
peak in P((Ni)) and polaritonic bunching (the extreme 
peaks in P((Ni}). This suggests that the system behaves 
very similarly to an Anderson-Mott insulator [27] . 
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FIG. 8. Probability distributions induced by disorder in 
the matter- light coupling g. (Top) Renormalized histogram 
P((Ni)) of the site occupation number as a function of disor- 
der. (Middle) Renormalized histogram P{E) of the site-site 
entanglement as a function of disorder. (Bottom) Renormal- 
ized histogram P((Z)) of the total atomic occupation number 
as a function of disorder. The light hoping is A = lg and de- 
tuning A = 5g. 

The presence of superfluid, Mott and Anderson-like 
states for intermediate disorder in the coupling g sug- 
gests that glassy phases would be induced in larger sys- 
tems of the Jaynes-Cumming-Hubbard type. In fact, a 
different situation was analyzed in [28 , in which there is 
disorder in the number of impurities per cavity. Since the 
number of atoms fluctuates then the intensity with which 
light couples to matter also fluctuates, and it was shown 
in [28] that such disorder induces glassy phases. Interest- 
ingly, even small versions of the quantum system present 
evidence of many diverse phases and behaviors expected 



only for larger quantum systems, which we discuss in the 
next section. 



IV. DISORDERED LARGE QUANTUM 
SYSTEMS 



To address the physics of large disordered quantum 
systems we apply a recently developed technique, namely 
stochastic- mean-field-theory (SMFT) [20]. This method 
has been shown to provide appropriate descriptions for 
the effects of disorder without over estimating coherence 
and fluidity and has already been successfully applied to 
the disordered Bose-Hubbard model [20]. The main rea- 
son for the effectiveness of the method is self-consistently 
determining the probability distribution for the mean- 
filed parameter P(a) (instead of a itself) through an it- 
erative process. 



A. Stochastic-Mean-Field-Theory 

Firstly we describe how to account for any on site dis- 
order (with constant photon hopping A), afterwards we 
describe the special case of hopping disorder. In the 
mean-field description every site has a number z of near- 
est neighbors. The mean-filed hamiltonian for the kth 
site depends only on the scaled sum 



Vk 



{k,j) a i- 



(7) 



The probability distribution for r] (we drop the site 
index for convenience) can be found from a simple and 
fundamental relation known as the convolution theorem 



Q(v) 



J ' ■■■ / II da i P ( a i) S ( V - A a 3 



(8) 



which reduces to the Fourier transforms (p(f3) = 
/ daP(a)e i(3a and 



-ir](3/A 



(9) 



The first step in the algorithm is to choose a trial dis- 
tribution for a (different from a delta centered at a = 0) 
and the desired distribution (in our case a Gaussian) for 
the disordered parameter £ (the detuning or matter-light 
coupling). Then we assume all ctj to be independent form 
each other and we determine the self-consistent distribu- 
tion 

P(a) = J J 'dq(ri)dp(Z)5(a-(a)), (10) 

such that (a) = r/]|a|G[£, 77]), with dq(rf) = drjQ(rj). 
The procedure is iterated until we observe convergence, 
that is, until pW (q,) i n the ith step is statistically close to 
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P( z+1 )(a). Finally, the average state of the site is given 
by the disorder- induced ensemble 



dq(v)dp(Z)\G[Z,Tj\){G[Z,rj\\. 



(11) 



To account for disorder in the photon hopping we must 
add another step to the procedure. It is convenient to 
work with the variable <\> = Aa whose probability distri- 
bution is given by (we use a subindex to distinguish the 
various distributions) 



P O) = J J dAdaP A (A)P a (a)S(Aa 



4>Y (12) 



Then we can determine Q{rj) through the usual Fourier 
transforms <p(P) = J d^P^e^ 



Q(v) 



1 

2^ 



-irjj3 



(13) 



and the procedure follows as for the case of on site dis- 
order. 



V. DISORDER-INDUCED TRANSITIONS 

The results presented in this section will allow us to 
conclude that the asymptotic effects of disorder (very 
large disorder) in the thermodynamical limit are very 
similar to the effects in small samples of the system. How- 
ever, there are some significant quantitative differences, 
for instance, there are in fact phase transitions induced 
by disorder in the thermodynamical limit. It should also 
be pointed out that our approach is a bit different in this 
section. From now on we work with single-site mean- 
field-theory rather than cluster- mean-field-theory, and we 
follow this strategy to avoid higher computational de- 
mands. This limits the applicability of the method and 
quantities like the site-cluster entanglement are no longer 
addressable. Nonetheless, we are able to increase the lo- 
cal effective dimension of the oscillator to 20 states. An- 
other difference between the approaches for small and 
large samples is that in the first case we fix the number 
of excitations in the system and it does not change as 
disorder increases. This is not the case in the present 
section, and in fact, the total number of excitations may 
change as a function of disorder. 

Using the SMFT approach we are able to perform an 
analysis of the thermodynamical limit. The method al- 
lows us to recover the probability distributions (under 
the single-site-mean-field approximation) for the various 
quantities we analyze to characterize the system. The 
average mean field parameter, for instance, can be read- 
ily evaluated as (a) = J aP(a)da. We follow the same 
ordering of presentation of the results: Firstly, we show 
the results for the detuning disorder, then for the hopping 
disorder and finally for the matter-light coupling. Given 
the unlimited nature of the disorder distribution we an- 
alyze (Gaussian) it follows that the insulating phases we 



present below are of glassy nature. Such phases have 
non vanishing number variance as opposed to the Mott- 
insulating phases [19] . Therefore glassy insulators can be 
characterized by vanishing superfluidity and non vanish- 
ing compressibility (which can be related to the number 
variance). However, we do not show the compressibility 
of the system, since the result can be readily anticipated. 
The compressibility increases in the insulating regions, 
thus glassy phases are established. 



1 2 3 



0.2 0.4 0.6 0.8 




10" 10" 




10 10" 



FIG. 9. Stochastic mean field results for disorder in the 
matter-light detuning with 5(A) = O.lg and (A) = 0. (Left) 
The average mean field parameter. (Right) The entanglement 
of the average state. 
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FIG. 10. Probability distributions as a function of disorder in 
the matter-light detuning with A — 10 -1 ' 9 g, (n — uj c )/g — — 1 
and (A) = 0. Black corresponds to vanishing probability. 
(Top-Left) Distribution for the mean field parameter, (Top- 
Right) for the entanglement of the average state, (Bottom- 
Left) for the cavity excitation and (Bottom-Right) for the 
atomic excitation. 
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A. Disorder in the matter- light detuning 

As show in figure ([9| the net effect of disorder in the 
detuning is to induce insulating behavior, indicated by 
the destruction of the fluid phase surrounding the Mott- 
lobes, in fact, the lobe structure disappears for significant 
amounts of disorder. The in-site entanglement shows 
that the system remains in superpositions of light and 
matter excitations for intermediate values of disorder. 
However, as we increase disorder the in-site components 
are either highly entangled or unentangled with higher 
probability. We can see the distributions as functions of 
disorder in figure (10). The transition from fluid to in- 



sulating is evident in the distribution of the mean field 
parameter. All this corroborates the small system pre- 
dictions and once again the system is mainly photonic 
for strong disorder. Interestingly, in the present limit the 
distribution of cavity excitation (P(N) in figure (fToj) ) is 
a series of delta functions (with different weights) cen- 
tered at integer values of the mean occupation, which is 
in agreement with the insulating and bunched behavior. 




0.2 0.4 0.6 




FIG. 11. Stochastic mean field results for disorder in the hop- 
ping with 6(A) = O.lg and A = 0. (Left) The average mean 
field parameter. (Right) The entanglement of the average 
state. 



B. Disorder in the photon hopping 

The effects that hopping disorder produces in the sys- 
tem are opposite to those produced by detuning disorder, 
as it is the case for small systems. In the current case 
disorder induces induces a fluid phase, as we can see in 
figure (11), and decreases the in-site entanglement indi- 
cating that the system becomes more photonic in nature. 
As in other situations analyzed, by looking at the distri- 
butions of the physical quantities as a function of disorder 
(not shown), we were able to observe the transition from 
insulating to fluid behavior as disorder in the hopping 
parameter increased. 
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FIG. 12. Probability distributions as a function of disorder 
in the matter-light coupling with A = 10 -1-9 g, (/i — oo c )/g — 
— 1 and A = 0. Black corresponds to vanishing probability. 
(Top-Left) Distribution for the mean field parameter, (Top- 
Right) for the entanglement of the average state, (Bottom- 
Left) for the cavity excitation and (Bottom-Right) for the 
atomic excitation. 



C. Disorder in the matter-light coupling 

Finally, the effects of disorder in the light-matter cou- 
pling, once again, resemble the ones induced by the de- 
tuning disorder. Even though both the detuning and 
coupling disorders induce insulating behavior, it should 
be pointed out they do it through very different physical 
mechanisms. 

In the detuning case cavities may be in lower frequen- 
cies in many sites which allows for more photons to local- 
ize (even several photons per cavity) . In the coupling case 
the strength of the polaritonic nature may be increased 
in cavities that are strongly coupled to their correspond- 
ing matter components inducing Mott behavior. And 
as we can see in figure ([l2] bottom left) in comparison 
with fig. ( 10 bottom left) the distribution for the cavity 
population does not present the higher order peaks, only 
the ones corresponding to zero or one polariton per cav- 
ity. This behavior is due to the fact that the sites with 
strong matter-light coupling prevent the accumulation of 
larger numbers of particles per site (Mott mechanism). 

Adding the information provided by the entanglement 
and atomic population distributions, we have found that 
a fraction of the sites assumes the Mott behavior and 
the rest are localized or even empty. Thus, we corrob- 
orate the small system analysis that suggests that the 
system behaves very similarly to the Anderson-Mott in- 
sulator [27]. It is, however, strikingly interesting that in 
one case (disorder in the detuning) it is the low frequency 
sites that dominate the resulting statistical behavior and 
in the other case (disorder in the matter-light coupling) it 
is the strongly coupled sites that dominate, even though 
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the distribution of the disorder parameter is gaussian and 
unbiased. 



VI. CONCLUSIONS 

We were able to characterize the phase diagram and 
the Mott-Superfluid transition of small and large sam- 
ples of the Jaynes-Cummings-Hubbard Hamiltonian us- 
ing entanglement measures between the various possible 
partitions of the components of the system. In partic- 
ular, we showed that these non-local measures identify 
more clearly where the transition happens. Furthermore, 
and more importantly, we also showed that entangle- 
ment measures distinguish which type of excitation dom- 
inates each phase which in turn allowed us to identify a 
crossover that is particular to this Hybrid system and 
does not have a purely bosonic analog. This behavior 
splits the superfruid regime into two, the first one domi- 
nated by polaritons and the second is purely photonic. 

For the disordered system we have shown that the sim- 
ple statistical treatment of small systems can be quite in- 
structive and allows us to draw conclusions that can be 
corroborated by the large system limit. We have shown 
that disorder both in the light-matter detuning and light- 
matter coupling induce insulating phases, however they 
do it through very different physical mechanisms. The 
former allows for photonic localization and bunching, the 



latter induces Mott behavior in a fraction of the sites that 
prevents the bunching. Furthermore, the cavity-cavity 
coupling disorder induces a glassy fluid phase. The rich 
in-site structure of the system leads to these diverse dis- 
ordered phases with very different statistics and physical 
meanings. 

A great amount of work remains to be done on the 
characterization of the JCH system. As a valuable point 
we suggest that an appropriate and efficient method 
should be applied to the study of the large hopping limit 
(with and without disorder) in which the mean field ap- 
proach adopted here is limited by the truncation of the 
state space. 

Finally, it is worth mentioning once again the meso- 
scopic aspect of the systems proposed to implement the 
JCH Hamiltonian as well as the increasing ability to ma- 
nipulate the different parameters of these systems, some- 
times even at an individual level. These properties sug- 
gest that it will be possible to carry on a thorough ex- 
perimental investigation of the effects of disorder and its 
relation to phase transitions and entanglement in many 
body physics in the near future. 
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